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We study the perturbative evolution of the static configurations, quasinormal 
modes and quasi normal ringing in the Apostolatos - Thorne cylindrical shell model. 
We consider first an expansion in harmonic modes and show that it provides a com- 
plete solution for the characteristic value problem for the finite perturbations of a 
static configuration. As a consequence of this completeness we obtain a proof of the 
stability of static solutions under this type of perturbations. The explicit expres- 
sion for the mode expansion are then used to obtain numerical values for some of 
the quasi normal mode complex frequencies. Some examples involving the numeri- 
cal evaluation of the integral mode expansions are described and analyzed, and the 
quasi normal ringing displayed by the solutions is found to be in agreement with 
quasi normal modes found previously. Going back to the full relativistic equations 
of motion we find their general linear form by expanding to first order about a static 
solution. We then show that the resulting set of coupled ordinary and partial dif- 
ferential equations for the dynamical variables of the system can be used to set an 
initial plus boundary values problem, and prove that there is an associated positive 
definite constant of the motion that puts absolute bounds on the dynamic variables 
of the system, establishing the stability of the motion of the shell under arbitrary, 
finite perturbations. We also show that the problem can be solved numerically, and 
provide some explicit examples that display the complete agreement between the 
purely numerical evolution and that obtained using the mode expansion, in partic- 
ular regarding the quasi normal ringing that results in the evolution of the system. 
We also discuss the relation of the present work to some recent results on the same 
model that have appeared in the literature. 

PACS numbers: 04.20.Jb,04.40.Dg 



I. INTRODUCTION 



The general analysis of the evolution of material systems interacting with their own 
gravitational field is undoubtedly a very important and very difficult problem in General 
Relativity. As usual in any physical theory, one can gain insight into the more difficult 
issues by considering idealized simplified models that, nevertheless, still preserve some of 
the relevant features of the general problem. One of these simplified examples was described 
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some time ago by Apostolatos and Thorne [1] . The material contents in this example is given 
by a cylindrical shell of counter rotating particles, and one is interested in the dynamics 
that results as a consequence of its interaction with its own gravitational field. As shown 
in [lj], once some appropriate choice of coordinates is made, it is a straightforward matter 
to obtain a set of coupled equations for the evolution of the matter and field variables, 
which is described in some detail in the next Section. In their original paper Apostolatos 
and Thorne were interested in general features of the evolution, in particular the possibility 
of naked singularities forming as the result of the collapse of these structures, but the 
detailed evolution of the shell in a situation close to one of its static configuration was only 
qualitatively mentioned, with no detailed results. This problem is interesting in particular 
because in the Newtonian approximation the shell is either in a static configuration or 
performs periodic oscillations about this static configuration, but, once the radiative modes 
of the gravitational field are introduced, one expects these oscillations to be eventually 
damped as the results of the transfer of energy to the radiative modes. Thus, in the fully 
relativistic description one expects to find the characteristic "Quasi Normal Ringing" (QNR), 
related to the "Quasi Normal Modes" (QNM)@ of the system. 

A recent analysis of the dynamics of the Apostolatos and Thorne model was given by 
Hamity, et.al. [3|. We notice, however, that the approximations introduced in the explicit 
examples considered in (3} are such that the radiative modes are effectively neglected, and, 
therefore, the results obtained are possibly relevant only close to the Newtonian limit. A 
detailed analysis of the full relativistic equations in the linearized approximation was carried 
out by the present authors in Q, where we obtained the general solution for the harmonic 
modes with real frequencies, and considered the solutions obtained by_general linear com- 
binations of these modes. Although in all the examples considered in |4( we found a stable 
evolution under perturbations of (essentially) compact support, at least two questions re- 
mained unanswered. The first was the possible completeness of the mode expansion, and 
its relation to the initial value problem for the system. The second was the apparent lack 
of QNR in the example solutions. The main purpose of the present work is to give answers 
to those questions. We first review in Sections 2 and 3, mostly for completeness, the main 
features of the model. Then, in Section 4, we review the derivation of the periodic solutions 
in the linearized regime, introducing a modified form, as compared with that given in Q, 
that allows us to establish, in Section 5, the completeness of expansion in the context of the 
characteristic value problem for the system, and, therefore, the completeness of the mode 
expansion. This completeness, in turn, as discussed in Section 6, implies the linear stability 
of static configurations under this type of perturbations, a result that differs from some con- 
clusions reached in jljj], and recent results obtained by Kurita and Nakao (6[. Next, in Section 
7, we discuss the existence of QNMs for the system. As shown there, these are related to 
the zeros of a complex expression that appears in the relation between the incoming and 
outgoing wave amplitudes as functions of the frequency. This expression has a complicated 
form involving Bessel functions that makes it very difficult its analysis. For this reason we 
have not been able to obtain explicit expressions for the frequencies of the QNM , or even to 
ascertain whether there is finite number or infinite number of these frequencies. Neverthe- 
less, by resorting to numerical methods we found at least a couple of QNM corresponding 
to a range of parameters characterizing static configurations. The values obtained are also 
displayed in Section 7. In Section 8 we consider the numerical evaluation of the integrals 
obtained in Section 7 that provide the evolution of the system for arbitrary characteristic 
data. The examples chosen display clearly the quasi normal ringing associated to one of the 
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QNM. The reason why this QNR was not seen in [4| is discussed in Section 12, using the 
Price-Husain toy model as an illustration. 

In Section 9 we consider again the full relativistic equations and assuming a general linear 
departure from a static configuration obtain the general linearized equations for the system. 
Then, in Section 10 we prove the existence of an associated positive definite constant of 
the motion that puts absolute bounds on the dynamic variables of the system, establishing 
the stability of the motion of the shell under arbitrary, finite perturbations. In Section 11 
we show that the corresponding set of coupled ordinary and partial differential equations 
for the relevant dynamical variables can be solved numerically as an initial plus boundary 
values problem, without resorting to expansions in terms of Bessel or other functions. We 
provide a couple of examples, using the same shell parameters as in the mode expansion, 
but a different form for the incoming wave form. The results, that clearly display quasi 
normal ringing are in complete agreement with those of the mode expansion, and with the 
computed QNM frequencies, although the techniques used in both approaches are totally 
different, confirming their separate correctness. 

We end the paper with some comments and conclusions, as well as comments on related 
work by other authors, in particular to that of reference 



II. THE APOSTOLATOS - THORNE MODEL 

The Apostolatos - Thorne model [l| describes the dynamics of a self gravitating cylindrical 
shell of counter rotating particles. Both the inner (M~) and outer (M + ) regions of the shell 
are vacuum space times with a common boundary E. The corresponding metrics may be 
written in the form, 

ds 2 ± = e 2 ^~ 2 ^ (dr 2 - dt 2 ± ) + e 2 ^dz 2 + e" 2 ^r 2 # 2 (1) 

where the (+) sign corresponds to the outside and (— ) to the inner regions. The functions 
ifj, and 7 depend only on r, t and satisfy the equations: 

VV + -i>,r ~ ^,tt = (2) 

Tjt = 2r^,i , 7,r = r [(</v) 2 + (^) 2 ] (3) 

The shell is located on the hypersurface E given by r = R(t), where r is the proper 
time of an observer at rest on the shell. We may interpret ip(r,t) as playing the role of 
a gravitational field whose static part is the analogue of the Newtonian potential. The time 
dependent solutions of (|2J) represent gravitational waves (Einstein- Rosen). Equation (J2J) is 
the integrability condition of Eqs. Q. The coordinates (z,<j),r) and the metric function 
ip are continuous across the shell E, while t and the metric function 7 are discontinuous. 
Smoothness of the spacetime geometry on the axis r = requires that 7 = 0,^ finite at 
r = 0, and dip/dr\ r=0 = 0. The junction conditions of M~ and M + through E require the 
continuity of the metric and specify the jump of the extrinsic curvature compatible with 
the stress energy tensor on the shell. The induced metric on E is given by 

dsl = -dr 2 + e 2 ^dz 2 + e- 2 ^R 2 d<f) 2 (4) 

Here V's(t) = ip+(R( T ), ^+( r )) — 4 , -{R( T )^-( T ))- The evolution of the shell is character- 
ized by R(t), which is the radial coordinate r at the shell's location and r the proper time 
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of an observer at rest on S. If we assume that the shells is made up of equal mass counter 
rotating particles, the Einstein field equations on the shell may be put in the form, 

2A 

i>%-^- 7 == m (5) 



4A\/i? 2 + e 2 ^ P 

X + -X' = - V „ (6) 

R 

where the constants A and J are, respectively, the proper mass per unit Killing length of the 
cylinder and the angular momentum per unit mass of the particles. The other quantities in 
( )5|6|) are given by, 

dt- 



X ± = ^ = +Ve- 2 (T±-^) + R 2 (7) 
or 

x(,% = i,±X± + iP±R (8) 
where a dot indicates a r derivative, and we also have, 

R 2 il>-X~ \R 2 X" J 2 e 2 ^X-X+ 

+ R 2 + e 2 ^J 2 ~ (R 2 + j2)3/2 + R ( R 2 + e 2fe J2) ( 9 ) 

These equations together with (121131) determine the evolution of the shell and of the grav- 
itational field to which it is coupled. 



III. STATIC SOLUTIONS 

We briefly review the conditions for a static solutions, corresponding to a shell of constant 
radius R{t) = Rq. We restrict to the case of a regular interior. In this case the interior 
must be flat and we may take 7_ = 0, ip- = 0, implying ip-z = 0. For the exterior field the 
general static solution (satisfying ip^, = 0) is of the form, 

^+(r) = -K\n(r/R ) , 7 +(r) = 7o + k 2 In (r/R ) (10) 

Then, 

X + (r) = e" 70 , X~ = 1 (11) 
Since R = 0, ipY, = 0, and ip~ n = 0, we find, 

k = 2^ , 1 q = 2\u[(R 2 + 2J 2 )/R 2 ] (12) 

and, 

_ j'Rq^pTr 2 

X ~ (2J 2 + R 2 ) 2 (13) 

It is important to notice that this relation can be satisfied for real Rq, J and A only if A < A c , 
with 

a/3 + Vl7 (Vl7 - 1) , s 

= v v yv ,_ = 15g79 14 

4(9 + ^) 
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Therefore, the system admits static solutions only if A < A c . 

It will be convenient to write J = xRq. Replacing in (fl3|) . we find, 



xWTTx^ , . 

X - (2x* + l)^ (15) 

which we now consider as an equation for function of A. 

It is now easy to check that there are no real solutions for x if A > A c . For A = A c 
we have a double solution x = x r = (vVl7- l)/2 = 0.88361 .... The interesting region 
for our analysis is the range < A < A c where we have two real and positive solutions for 
x, one larger and one smaller than x c , which approach respectively +oo and as A — > 0. 
This implies that for appropriate A and J we have two static solutions. It was concluded in 
previous analysis jl] that, at least from a perturbative point of view, one of these solutions 
(the one with the larger R ) would be stable and the other (with the smaller Rq) would be 
unstable. This conclusion, as shown in the following Sections, is not supported by the present 
analysis of the perturbative evolution that follows from initial data in the neighborhood of 
either static configuration. In the next Section we derive the appropriate equations for a 
perturbative analysis, and then show how this construction can be used to solve the evolution 
equations given characteristic data. 



IV. LINEARIZED PERIODIC SOLUTIONS WITH A REGULAR INTERIOR 

In this Section we consider the problem of finding linearized periodic solutions for our 
system imposing the condition of a regular interior. This problem was already considered 
in where the main idea was to obtain stationary periodic solutions, but left open the 
problem of relating these solutions to the initial value problem. Here we will consider a 
slightly modified approach, where we express the solutions in M + in terms of incoming and 
outgoing wave solutions for the radiative part of 

We first notice that if we impose regularity on the symmetry axis r = 0, and assume 
that the interior region is empty but may contain gravitational radiation, then %p_ admits 
periodic solutions of the form |5j, 

^_(t_,r) = AiJ (|fi-|r)e in - t - (16) 

where Jo is a Bessel function, and we consider A± as a quantity of first order in perturbation 
theory. Then, restricting again to linearized order we may set, 

T _(t_,r) = (17) 

since, in this case, from ([3]), 7_ is of second order, and, therefore, also to the appropriate 
order, we may also set, 

t-(r)=r (18) 

We assume a periodic perturbation around an equilibrium configuration characterized by 
Ro and J and therefore take, 

R(r) = R + ^e mT (19) 
where we assume that £o is a l so of first order. 



6 



We need to specify now a corresponding solution for ip + . In our previous analysis, [4|], we 
noticed that the general periodic perturbation for ip + may be written as, 

MU,r) = -/tin + [A 2 J (\n + \r)+B 2 Y (\n + \r)]e in ^ (20) 

where Y is a Bessel function and A 2 and B 2 are first order quantities. To this order we 
have, 

/ f \ 

.2i„ i o,. r /i 7 /in i„\ i d v /in i„\i „is2+*-i 



7+ (t+, r) = to + « 2 In ( — ) - 2k [A 2 J (\n+\r) + B 2 Y (\Q + \r)} e^+ 4 + (21) 



where k and 70 are given by ( fl2l . A straightforward calculation [4| shows that consistency 
of the equations at first order requires ft_ = ft, 



(R 2 + 2J 
? 4 



2\2 



ft+ = v — ^ft = e 7o ft (22) 

Rr, 



and, 



<+ = MrW T = e ~ 7 ° r (23) 

Replacing now in (jHJ), (EJ), and ([9]), and expanding to first order, we find a set of three linear 
independent equations for A±, A 2 , B 2 , and £ 0; f° r every choice of ft, Rq and J, and, therefore, 
three of these quantities can be solved in terms of the fourth. Since all the evolution equations 
are linear, we may consider linear superpositions of these solutions to obtain more general 
(non periodic) solution. In Q, for simplicity, we chose £0 as the independent variable, and 
obtained corresponding expressions for the other three amplitudes. Assuming an arbitrary 
dependence of £0 on ft, given by the function £(ft), R(t) is given by, 

/+00 
£(ft) exp(iftr)dft (24) 
-00 

while for ip_(r,t_) and tp + (r,t + ) we obtain expression in the form of Fourier transforms of 
£(ft) multiplied by complicated expressions involving Bessel functions of rQ± and Rofl±. 
Details are given in Q. We notice that although in principle, (I24p is completely general, 
and one has a definite £(ft) for any possible R{t), in practice, what one would like to solve is 
the problem of the evolution of the system given some initial condition, and because of the 
structure of the equations that resulted in [ij , the best one could do was to assume some form 
for £(ft), and find the corresponding evolution of the system. Carrying out this program it 
was found there that there are rather general evolutions, in the sense that they represent 
the interaction of the shell with an incoming gravitational pulse of rather arbitrary shape, 
for which the response of the shell is given precisely by £(ft). An important question here is 
actually how general these solutions of the problem are, namely, whether they include every 
possible evolution of the system, or on the contrary, they only represent a restricted set. 
Since we had no proof of the completeness of the set of modes, this question is highly non 
trivial. Further, the scheme developed in Q did no appear to show the quasi normal modes 
and quasi normal ringing expected in a system such as this, that has a simple Newtonian 
limit, where the shell can execute periodic motions. One would expect that the presence 
of radiative modes in the general relativistic version should lead to damped oscillations, at 
least in some appropriate limit, and this was not immediately apparent in the discussion 
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given in |4|]. As we shall show here, both the completeness problem and the question of the 
QNR can be answered by a simple change in the formalism presented in 0]. This can be 
achieved as follows. The expression, 



ip in (t+,r) 



F e i: 



M\Sh\r) 



10, 



■Y (Wr) 



(25) 



where F and Q 2 are constants, is a solution of (j2J) for For large r, using the asymptotic 
expansions for the Bessel functions J and Y , we find, 



ipin(t+,r) 



1- &\ 

ito 



(26) 



and, therefore, ( |25j) represents a purely incoming wave solution of fl2J). 
Similarly, 

Vw(t+, r) = G e m *+ J (\n 2 \r) - i^Y (\n 2 \r) 

"2 



(27) 



where G is a constant, is a solution of (j2j), for -0 + , which, for large r has the asymptotic 
behaviour, 



ipout(t+,r) ~ G 1 



1 



7rr|f2g 



1 + 

1 60 



(28) 



and, therefore, ( 127)) represents a purely outgoing wave solution of (|2j). We notice now that 
we have, 



(F + G)J Q W 2 \r) + i—^r(F - G)Y (\n 2 \r) 

it2 



^ n (t + ,r) + iJ out (t + ,r) 
Therefore, instead of ( 120]) we may write the full solution for ip + (t + ,r) in the form, 



(29) 



if)+(t +i r) = -/tin — +^ in (t + ,r) + ^ out (t+,r) 



where we have the identifications: 



A 2 = F + G 



(30) 



5o 



2 



(F - G) 



(31) 



In what follows we will consider F and G as first order quantities. To this order, from 
), we have, 



7+(t+,r) = 70 + K 2 ln ( — J - 2k (ip in (t + , r) + ip out (t + , r)) 

H 



(32) 



where «; and 70 are given by ( 1121) . The important point here is that we have now expressions 
for tp + and 7+ where we have separate amplitudes for the incoming and the outgoing wave 
parts. 
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Notice that since ( 130]) is just a reordering of terms in (1201) . all the previous relations 
regarding t±, Q, and Q± are valid. Then replacing again the expressions for the dynamic 
variables, but using now (130]) and ( 132|) . in flSJ), (]6]), and fl5J), and expanding to first order, we 
find this time a set of three linear independent equations for A±, F, G, and £o- I n this paper 
we will be mainly interested in the evolution of an initially static shell, under the influence 
of a bounded incoming compact pulse. We choose therefore to solve these equations for Ai, 
G, and £o i n terms of F. The resulting explicit expressions are, unfortunately, rather long 
and difficult to read. Nevertheless, they have in general the structure, 

* = §F 

G = §F (33) 

6, = ^ 

where Hi , -i = 1,2,3 are bounded regular functions of (Ro, J, Q) and D is given by, 

D = - [2 [R 2 + J 2 ) J 2 R ((4 J 4 + QR 2 J 2 + R 4 ) Vt 2 - 2 J 2 ) J R ) 

- (i? 2 (2 J 2 + Ro 2 ) 2 {Ro 2 + J 2 ) 2 ti 2 - J 2 (2 i? 6 + J 2 (2 J 2 + Ro 2 ) 2 )) OJi (|fi| ^o)~ 

J (|J) 2 |i2o)- ^iod^l^o)) 

- ~(V (i? 2 + J 2 ) 2 fi 2 - J 2 (2R 2 + J 2 )) Jo (M R ) 

+2J 2 nR (R 2 + J 2 ) Ji JSo)] (2 J 2 + i? 2 ) 2 

(ft Ji (|fi 2 | i? ) - ^ n (|n 2 | iJo)) (34) 
where Q 2 — ft+- 

In spite of its appearance it is not difficult to show that D has no zeros for real Q. A 
zero of D for real Q would imply that the imaginary part of 



\n\j (\n, 


> Ro) — 


inY (\n 2 


Ro) 




R )-i 


0| Yi (|fi 2 


Ro) 



(35) 



vanishes. But this is given by, 



j (\n 2 


R )Y 1 (\n 2 


\Ro)-Ji(\n 2 


Ro)Y (\n 2 


Ro) _ 




2 








Ro) 2 


tiRq 


n 2 \(Ji(\n 2 


Ro^ + Y^ 


Ro) 2 ) 



(36) 

where we have used properties of the Bessel functions to obtain the right hand side, and, 
therefore, the imaginary part is non vanishing for all Q. 

This result is crucial, because it means that if F depends regularly on O, the solutions 
obtained using ( 133]) and ( 134]) are well defined modes for all Q in —00 < Q < +00. Therefore, 
because of linearity, arbitrary linear combinations of these modes will also be solutions of 
the problem. Before discussing the possible completeness of this set of modes we will write 
explicit forms for the expressions we have in mind. These are given by, 



ipin(U,r) = / F(Q)e 



i'2 



dn (37) 
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where F(Q) is an arbitrary complex function of Q which we will assume at least square 
integrable.Then, taking into account ( )33|) . we have for the remaining dynamical variables, 



tp out (t + ,r) 
*(r) 



Hi 
D 



F(n)e im -j (\n\r)dn 



+00 H 7 



D 



F(D,)e 



oc 



"2 



dQ 



D 



F(Q)e inT dQ 



(38) 



We shall analyze the properties of this expansion in the next Section, and show that they 
solve a well defined characteristic value problem for our system, and, therefore, in this sense, 
the mode expansion is also complete. We remark, for completeness, that (137|) . for real Q, 
represents the most general solution for ij; in (t + ,r), such that both ij; out (t + ,r) are bounded 
for large r. 



V. THE CHARACTERISTIC VALUE PROBLEM 



Consider again fl37j) . The problem that we have in mind is one where at some given 
time the shell and a large neighbourhood of the space time that includes the shell and the 
symmetry axis, are in a static state, but there is a finite gravitational pulse incoming from 
large r. In other words, if i?o corresponds to a static solution, then for some large time in 
the past and some R s » Rq, the solution in < r < R s coincides with the static solution. 
For r > R s , on the other hand, there is region where ipi n (t + ,r) is non vanishing. To make 
this more precise, we recall that for large r an incoming wave solution for ip + (t + ,r) has the 
asymptotic form, 

Vin(*+,r) ~ —f(t+ + r) (39) 



where f(u) is an arbitrary function of (essentially) compact support. We may compare this 
with the asymptotic form for large r of ()3j 



r 




L /(t+ + r)~/ F(n)J—^— ( 1 -i^t ) e in2it++r) dQ (40) 



which can be inverted to give, 

in 1 3/2 r+co 

F (Q) — ' 2 — / f(u)e- iQ2V du (41) 



Therefore, we find a one-to-one relation between F(Q) and characteristic data given for 
r — > oo, t + — > — oo. In this sense, we have proved the following result: The set of modes given 
by ( Tffffj) and (3$ ), used in the construction of (31) and (38\) is complete, and, in particular, 



(3l\ ) and §3E) provide the evolution of the corresponding dynamical variables for arbitrary 
characteristic data given at r — >■ oo ; t + — > — oo. 

In the next two Sections we shall use these results to analyze the stability of the static 
solutions and then consider the presence of quasi normal modes (QNM) and their associated 
quasi normal ringing (QNR). 
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VI. STABILITY OF THE STATIC SOLUTIONS 



Conceptually, we may say that a static configuration of the shell is stable if, given data 
arbitrarily close to that configuration, the system evolves towards that static configuration. 
Otherwise we would say that the static configuration is unstable. But, as we have already 
indicated, we have explicit forms for the evolution of the system when it is perturbed from 
an initially static configuration by an incoming gravitational pulse of (essentially) compact 
support. Each one of these pulses is uniquely characterized by a corresponding function 
F(Q). Concentrating in particular on £(r), we have, 



AitfRj (.J 2 (2J 2 + 3i? 2 ) 



7T 



R Z 



nR (2 J 2 

and, using ( B1J) , we can check that for large Q 



J Ro) - {R 2 + J 2 ) Ji(P\Ro 



(42) 





H 3 




2V2R 7/2 cos(fii? + vr/4) 


1 




D 




v/?(2J 2 + /V) (Ro 2 + J 2 ) 


|Q| 3 / 2 


plus terms of order \Q\ 


5/2 


As 


consequence, on account of 


(EHD, a 



(43) 



F(Ti), £(r) is the Fourier transform of a square integrable function for any J and Rq, but 
then, we must have £(r) — > for r — > +oo, and, therefore, all static configurations are stable. 
More explicitly, we have shown that all finite admissible gravitational pulses incoming from 
large r can be described by an appropriate function F(Q), and that this function is in one- 
to-one correspondence with the shape of the pulse at large r. We have also shown there is 
a unique evolution for the system corresponding to each function F(Q). This evolution is 
such that the dynamic variables are essentially given by the Fourier transform of a square 
integrable functions and, therefore, they are themselves square integrable and must vanish 
in the limit |r| —> oo. But this shows that all initially static configurations of the shell that 
are perturbed by an incoming pulse of the type described above, will eventually settle back 
to their original static configuration, and, in this sense, they are all stable. 

Just for completeness, it should be clear that we have not found a self adjoint extension 
associated to the coupled partial plus ordinary system of differential equations of this 
problem, so that we could not state completeness of the mode expansion in the same way 
as when such an extension is possible. What we have found is that the characteristic 
value problem can indeed be solved by our mode expansion, since it reduces essentially 
to a Fourier transform, and, that in this sense, the expansion is complete, and proves 
that the static configurations of the system are stable under a perturbation that has 
the form of a bounded, (essentially) compactly supported incoming pulse. But this is 
essentially all that is required to consider the system physically stable. It is still an open 
and interesting question to find the relation, if it exists, between the system of equations 
governing the evolution of the system and an associated self adjoint problem that would be 
more useful an the analysis of the initial value problem. Nevertheless, as show in Section 
10, for the general initial plus boundary value problem we have a conserved, positive defi- 
nite constant of the motion, and therefore we have stable evolution for any finite initial data. 

We have already seen that for a shell with given values of A and J there may be just 
one, two or no static solutions of the equations of motion. It has been suggested |3j that in 
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the case of two solutions only one is stable and the other is unstable. However, our present 
analysis does not show a qualitative difference in the behaviour of two static solutions under 
perturbations. There appears to be, nevertheless, a more subtle difference between the 
solutions when we analyze the evolution in more detail. This is related to the presence of 
quasi-normal modes, and is considered in the next Section. 



VII. QUASI NORMAL MODES 



In this Section we analyze the presence of quasi normal modes and the associated quasi 
normal ringing in our system. Very roughly, QNM and QNR appear in systems coupled to a 
field that admits a decomposition in periodic incoming and outgoing waves. (See e.g. ^ for 
more details, and further references). Generally QNM are associated to non trivial solutions 
where the incoming wave amplitude vanishes. This requires, generally, a complex value for 
the period, and therefore, in general, the corresponding outgoing and other amplitudes are 
unbounded in space and time, and, therefore, are not physical. In our problem, we have 
seen that for the periodic solutions, the dynamical variables other than F(Q) are given 
by expressions of the form HiF/D, and therefore, QNM, that is, nontrivial solutions with 
F — 0, can exist only for values of Q such that D — 0. We have already seen that D 
is non vanishing for real Q, but it may vanish for complex values of Q. These zeros of D 
introduce complex poles in the expressions for £ and the other field variables whose complex 
frequency does not depend on F although the amplitude of their possible contributions 
does depend on F. Since £, A\ and G are given by Fourier transforms of these expressions 
containing complex poles, we may get contributions dominated by these poles. These are 
the QNR amplitudes associated to the QNM. Therefore, to find the QNM we need to find 
the complex zeros of D. In principle this is the simple matter of solving the equation D = 
for Q. Unfortunately, if we consider the explicit expression for D given by (134)) . we notice 
that it contains Bessel functions of different types and arguments depending on f2, so that 
only a numerical computation of the zeros appears as feasible, but even in this case the 
accurate computation of these zeros is a complicated task. In this paper we only attempt a 
rather rough computation of the zeros closest to the real O axis in the complex Q plane, as 
these would correspond to the most noticeable QNR of the system. 

For this explicit computation we go back to fl33|) and notice that if define x = J/Ro and 
E = QRq, the equation D = may be written in the form, 

= [2 (l + x 2 )x 2 (E 2 (4x 4 + 6x 2 + l) - 2a; 2 ) J (|E|) 

-E (S 2 (2x 2 + l) 2 (1 + x 2 ) 2 - x 2 (4x 4 + 2 + x 2 + 4x 6 )) J, (|E|)~ 

(44) 



(-|£| J (|E| (2x 2 + l) 2 ) +*EV (|£| (2x 2 + l)' 
- (EJi (JE| (2a; 2 + l) 2 ) - i |E| Y x (jE| (2a; 2 + l) 



E 2 (1 + x 2 ) 2 - x 2 (x 2 + 2) ) J (|E|) + 2Ea; 2 (l + x 2 ) J x (|E|) E (2a; 2 + l) 



Actually, since we are looking for zeroes near the real Q axis, we should make the re- 
placements |E| = E to look for zeros in the region 9?(f2) > 0, and |E| = — E in the region 
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K(fi) < 0. In any case, the zeroes of D satisfy the scale law, 

n QNM = -j-g (-£) (45) 

where Q is some complex function of its real argument. 

We searched for zeroes of D using two procedures. In the first we noticed that if we look 
for zeroes in the 9ft (E) > region we may replace |E| = E in (jUJ), and after canceling some 
common factor of E, we find a maximum power E 3 in the coefficients of the Bessel functions. 
We formally solve for this factor and obtain, 

E 3 = [(2 (l+x 2 ) ((4x 4 + 6x 2 + 1) E 2 -2x 2 ) J (S) - Ji (E) E (4x 4 + 2 + x 2 + 4x 6 )) 

(J (S (2x 2 + l) 2 )-zF (s (2x 2 + l) 2 )) 

- (2 J x (E) E (1 + x 2 ) - J (E) (x 2 + 2)) 

(iYi(E (2x 2 + l) 2 )-J!(s (2x 2 + l) 2 ))E (2x 2 + l) 2 ]x 2 (46) 
"((^i(s (2x 2 + 1) 2 )-J!(e (2x 2 + 1) 2 )) J (E) 

+ (Jo(e (2x 2 + 1) 2 ) -.T (e (2x 2 + 1) 2 )) J x (E)) (2x 2 + l) 2 (l + x 2 ) 2 ] ~* . 

We can use this equation in an iterative scheme where we input a value for E on 
the right hand side, and the cubic root of the (complex) number obtained is inserted 
again as an "improved" value for E. We have numerically checked that this works very 
well for sufficiently small values of x, say x <~ 0.5, but stops converging for larger 
values of x. A different, in a way more direct, although not very accurate method 
that we have used is to consider again (I44p . fix a value of x, and plot the real and the 
imaginary parts of (144 p . as functions of 9ft(E), for fixed 3(E). With this method one 
can easily visualize the possible common zeros of the curves and adjust, by trial an 
error, the best values of 9ft(E) and 3(E). We have found by this method that there are 
complex zeros of D at least in the range < x ~ 2, and possibly for larger values of x. 
We have found zeros both near 9ft(E) > and 9ft(E) < 0. Some of these are shown in Table I. 

The first thing to notice is that all these zeroes have positive imaginary part, as one 
would expect from the stability arguments of the previous Section. We also notice that 
there are zeros both for 3ft(E) > and !R(E) < 0, although having different imaginary parts 
and different |3ft(E)|. In general the imaginary part of E increases with x, and the modes 
become very strongly damped for x of the order or larger than one. All these results are 
in good agreement with the "first class solutions" of |6[ for x < 2, although they obtain 
a much larger set of solutions. Regarding this point we remark that, unfortunately, the 
analytic structure in the complex Q plane of the integrands in (1381) is very complicated, and 
the exact relation between these zeroes and the evolution of the dynamical variables is not 
easily established. To obtain more information on the subject of QNR for our system we 
carried out numerical integrations of £(t), as given by (1381) . for simple incoming F, and 
found examples of QNR, with parameters close to those of the zeros of D found here. This 
is detailed in the next Section. Next, in Section IX, we consider again the perturbation 
problem in general, and obtain a set of coupled equations directly for £, ^_ and without 
resorting to a mode expansion. This allows us to set up an initial value problem that can 
be solved fully numerically, without the use of the mode expansion or Bessel functions. 
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X 


K(S) > 


K(S) < 


0.10 


0.1396+ 0.00004i 


- 0.1408+ 0.00015i 


0.20 


0.2685+ O.OOlOli 


- 0.2760+ 0.0035i 


0.30 


0.3772+ 0.00575i 


- 0.3939+ 0.01775i 


0.40 


0.4612+ 0.01722i 


- 0.4840+ 0.0479i 


0.50 


0.5202+ 0.0366i 


- 0.5410+ 0.09M 


0.60 


0.5570+ 0.0628i 


- 0.5676+ 0.13998i 


0.70 


0.5750+ 0.0938i 


- 0.5704+ 0.1887i 


0.80 


0.5780+ 0.127i 


- 0.5552+ 0.2337i 


0.90 


0.5704+ 0.16125i 


- 0.5305+ 0.273i 


1.00 


0.5535+ 0.1943i 


- 0.4970+ 0.307i 


1.20 


0.5010+ 0.2542i 


- 0.4200+ 0.3588i 


1.40 


0.4350+ 0.3033i 


- 0.3330+ 0.395i 


1.60 


0.3590+ 0.342i 


-0.2350+ 0.4201i 


1.80, 


0.2736+ 0.37175i 


- 0.0960+ 0.4378i 


2.00 


0.1708+ 0.3947i 





TABLE I: The zeros of D. 



VIII. NUMERICAL EVALUATION OF THE INTEGRALS IN (l38|) 

We have carried out explicit numerical integrations of (1551) to obtain £(r), with the choice, 



1 3/2 



F ^ = (^y exp( - t!/6) (47) 



This corresponds to an asymptotically incoming pulse of the form, 

f(t + + r) = y/6e X p(~\t + + r\ 2 ) (48) 

As a first choice we took Rq = 0.25 and J = 0.125. This corresponds to x — 0.5 and, 
therefore, is within the range considered "stable" in previous analysis. The resulting £(t) 
is displayed in Fig. 1. We can see the characteristic shape of a QNR, that is, a damped 
oscillation. As usual, we can get a better look at this shape by displaying, as in Fig. 2, 
a = ln|£(r)| as a function of r (thick line). We also display in Fig. 2 an approximate fit 
using the function ln(| cos(ui r r + <p)\) — UiT + q, where u r = 2.08, <p = 0.01, = 0.16 and 
q = 4.28, although in the actual figure, for clarity, we set q = 2 to avoid superposition of 
the two graphs. These parameters are in good agreement with those corresponding to the 
QNM with SR(ft) > for x = 0.5 shown in Table I. 

As a second example we took R = 0.25 and J = 0.25. This corresponds to x — 1.0 and, 
therefore, is outside the range considered "stable" in previous analysis. The resulting £(r) is 
displayed in Fig. 3. We can see again the characteristic shape of a QNR, but, in this case, as 
a strongly damped oscillation. This result is also in good agreement with the fact that the 
QNM for x = 1.0 shown in Table I corresponds to a strongly damped oscillation. We remark 
once again that the analytic structure of the integrand in (155]) is far from simple, and that 
although several sets of complex QNM frequencies can be established rather accurately, their 
effect on the resulting amplitudes cannot be easily established. Nevertheless, the integrals 
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30 



-60 



FIG. 1: £(t) as a function of r for an incoming pulse of the form (|48p . for Rq = 0.25 and J = 0.125 

in (]38|) are taken along the real f2 axis, and therefore, they are insensitive to the manner in 
which the integrand might be extended to, e.g., the upper complex- Q plane, where the poles 
corresponding to the QNM reside. 

In the next Section we consider a different approach, where we linearize the equations 
of motion and obtain a set of coupled linear equations for the amplitudes, together with a 
set of boundary and matching conditions that make possible to establish (and numerically 
solve) a well defined initial data problem. 



We consider again the set of dynamical equations and matching conditions for our prob- 
lem, assuming that we are close to a static solution, characterized by certain values of J and 
R . We therefore write, 



IX. THE GENERAL PERTURBATION PROBLEM 



R(r) 



Ro + e £(t) 
e Xi(t_,r) 

-«ln(r/i2o) + e X 2 (t+,r) 



(49) 



with Xi satisfying the equations, 









(50) 



dt 2 _ 
d 2 X2 







(51) 



dtl 



dr 2 



r dr 
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FIG. 2: a = In |£(r)| for an incoming pulse of the form (|48p . for i?o = 0-25 and J = 0.125 (thick 
line). The thin line is an approximate fit by ln(| cos(2.08r + 0.01) | ) — 0.16r + 4.28. Notice that the 
actual plot of this fit has been displaced down by 2.5 to avoid cluttering the figure. 
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FIG. 3: ^(r) as a function of r for an incoming pulse of the form (|48p . for = 0.25 and J = 0.25, 
that is x = 1.0. Notice the strong damping of the oscillations after the arrival of the incoming 
pulse. 
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and expand the remaining equations to first order in e. Using these expansion we first obtain, 

7-(*-,r) = 0(e 2 ) (52) 
l+(t+,r) = 70 + K 2 \n(r/R ) - 2e k Xa (<+>»") 



In these equations k and 70 are given by ( I12j) and ( I13j) . so that the zeroth order terms 
coincide with the static solution of Section III. From these, again to first order we find, 



dt_ 

dr 
dt + 

dr 



1 + 0(e) 
e" 70 + 0(e) 



(53) 



The explicit form of the terms of order e is not required in what follows. Choosing appro- 
priately some integration constants we set, 



t_(r) = r + 0(e) 



(54) 



Ur) 



-70 



r + 0(e) 



Next we replace in the matching conditions and expanding again to first order we obtain 
two equations for Xi an d X2 on the boundary r = R 







2 J 4 (QR 2 J 2 + 4J A + 3R 4 ) 



(R 2 + J 2 ) (2J 2 + J R 2 ) 2 J R 



£(t) 



2i?n t/ 4 



(^o 2 + J 2 ) (2J 2 + J R 2 ) 2 
2 J 2 (4J 2 + R 2 ) 



(2J2 + R 2 yR 
2J 2 



X2 (e" 70 r 



Xi (t, Rq) 

Ro) 



dxi (r, r) 



<9/' 



(55) 



-Rn 



<9X2 (e" 70 r, r) 



(2J 2 + J R 2 )' 



<9r 



r=-R 



-Rn 3 



C{t)+ X i (r,R )-X2 (e-^r,Ro) 



and an equation of motion for £(r) 



dr 2 



J 2 (4 J 6 + 6R 2 J A + SJ 2 ^ 4 + 2i? 6 ) 



2 - 2 




(2 J 2 + J R 2 ) 2 (i? 2 + ^ 2 )^i 
J 2 ^ 3 (3J 2 + 2i? 2 ) 



{2.P + R, 2 ) z {R, 2 + .P) 
(4J 2 + Rq 2 ) R J 2 

(R 2 + J 2 ) (2J 2 + J R 2 ) 2 



Xi (t, -Rq) + 



^o 2 ^Xi (t, r) 



iV + </ 2 



<9r 



(56) 



r=-R 



X2 (e-^r,i2o) 



For explicit numerical integration it turns out to be convenient to introduce instead of 
X2(t+,r) a new function X3( r ? r ) defined by, 



Xa(T,r) = X2(e 7o r,r) 
and, therefore, satisfying the equation: 

(2J 2 + J R 2 ) 4 9 2 X3 (r,r) , <9 2 X3 (r,r) , l3 X3 (r,r) 



i2g 



«9r 2 



<9r s 



r <9r 



(57) 



(58) 
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This way all the dynamic variables involved evolve directly in terms of r. This simplifies 
the numerical treatment of the problem as it eliminates the need of separate grid spacing 
for r < R and r > R . But, most important, using % 3 we may demonstrate the existence 
of a crucial constant of the motion, as is shown in the next Section. 



X. A POSITIVE DEFINITE CONSTANT OF THE MOTION. 

As discussed in the previous Section, we have a complete solution for characteristic data, 
and we can use this fact to prove stability regarding the evolution resulting from that type 
of data. This, nevertheless, leaves still open the question of the evolution of general initial 
data. In other words, the possibility of the existence of initial data that somehow is not 
registered as characteristic data, but such that it renders the motion unstable, in the sense 
that, i.e. £(r) can acquire arbitrarily large values starting from finite initial data. To answer 
this question we consider again ( 1531) . replace \2 in terms of X3, and solve for xi{t,Ro) and 
X 3 (r,R ). We get, 



Xi ( r > #o 



and, 



X3 0, R Q ) 



R (Rq 2 + J 2 ) (2 J 2 



Ro 2 ) 2 



2.P (6 R 2 J 2 + 4 J 4 + RO 4 ) 
Rq 5 {Rq 2 + J 2 ) 



d , 



2.P (6Rq 2 J 2 + 4P + RO 4 ) \dr 
J 2 (4R0 2 J 2 + 4J 4 -R 4 ) 
R 3 (6R 2 J 2 + 4P + R 4 f (r) 



d , 



2P 



R (Rq 2 + J 2 ) (2J 2 + Rq 2 ) 
' 2 J 2 (6R 2 J 2 + 4 J 4 + RO 4 ) 

Rq 5 (Ro 2 + J 2 ) 
2 J 2 (6R 2 J 2 + 4J 4 + RO 4 ) 
J 2 (3R 2 + 2J 2 ) (2J 2 + R 2 ) 
Ro 3 (6R 2 J 2 + 4J 4 + Ro 4 ) 



Or 

d_ 

dr 



Xi {i~,r) 
X3 (r,r) 



Ur) 



r=R 



r=Ro 



r=R 



r=R 



(59) 



(60) 



Replacing now (1591 and (16"0"1) in (1561) . upon multiplication by d^/dr, and some rearrange- 
ment, we find, 



di{r)d 2 j{r) 2 J 2 dj{r) 

dr dP QR 2 P + 4J 4 + i? 4UrJ dr 

R 2 (4R 2 J 2 + 4J 4 -R 4 ) d^jr) f d_ \ 

2 (R 2 + J 2 ) (6 R 2 J 2 + 4 J 4 + .Ro 4 ) dr \dr Xl (T ' T) ) r=j 
(3 -Rq 2 + 2 J 2 ) R 6 c%(t) (d_ 

2 (R 2 + J 2 ) (6 R 2 J 2 + 4 J 4 + R 4 ) (2 J 2 + Rq 2 ) dr \dr 
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We notice that the first line (16T|) has the form of a total r-derivative. We may put the 
whole expression in this form as follows. We first take the derivative of (|59|) with respect 
to r, solve for d^/dr and replace it in the second line in (16T|) . Next, we derive (160]) with 
respect to r, solve again for d^/dr, and replace it in the third line in (|6"T|) . After a new 
rearrangement we get, 



d£d 2 £ 



2J 2 



drdr 2 6 Rn 2 J 2 + 4 J 4 + R dr 



(2 J 2 + R 2 ) 2 d 2 xi d 2 X A ({2J 2 + Ro 2 Y d X i d X 



2\ 2 



Rq 



ii 



drdr drdr 



dr dr 



4 (2 J 2 + i? 2 ) </ 4 (6 i? V 2 + 4 J 4 + i? 4 ) 

#Xi ^Xi ^o 9 dx 3 



(62) 



2 (i? 2 + J 2 ) J 2 dr Or 2 (R 2 + J 2 ) (2 J 2 + R 2 ) 2 J 2 dr dr 

where, again, all derivatives of Xi an d X3 are evaluated at r = Rq. 

The first two lines in ( 162]) are total r-derivatives. To analyze the last line we notice that, 
on account of ( |50l) we have, 



d 

dr 



Ro 



dxi 
dr 



dxi 
dr 



dr 



Ro 



dxi(r,r)dxi(r,r) 
dr dr 

dxi(r,R ) dxi(r,r) 



dr 



dr 



Ro 
o 

r=R 



(63) 



where in the last line we have assumed that xi satisfies the (regularity) boundary conditions 
that Xi( T ,Q) is finite, and that (dxi/dr)\ r=0 = 0. 
Similarly, from ( I58p . we have, 



d 

dr 



Ro 



(2J 2 



Rl 



dxs 
dr 



dxs 
dr 



dr 



9x3(r,r) dx3 (r,r) 



— Rq 



dr dr 
9x3(t,R ) 9x 3 (T,r 



Ro 



dr 



dr 



(64) 



r=R 



where we have assumed that ip 3 (r,r) is such that the integral exists, and, therefore, there is 
no contribution for r — > oo. 

Comparing ( 163]) and (164"]) with (162]) we find that, if ipi(r, r) corresponds to a solution 
of (158]) that is regular for r = 0, (and, therefore, r(dipi/dr)(dipi/dr) — > for r — )• 0), 
and ^(r, r) is such that r(dip3/dr)(dip3/dr) — > for r — >• oo, (which must happen for the 
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integral in ([6 



E H 



dr 



to exist), then the quantity, 

2 



+ 



P 



6 Rq 2 J 2 + 4 J 4 + Rc 



D 14 
tin 



(2J 2 + R 2 ) 2 d Xl d X3 



8 (2 J 2 + Rq 2 ) 2 J 4 {QR Q 2 J 2 + 4 J 4 + i? 4 ) 



i?0 



<9r <9r 



=-Ro 



2 (R 2 + J 2 ) J 2 



Ho 



<9r 



i?n 8 



2(i? 2 + J 2 ) (2J 2 + i? 2 ) i J 2 



flo 



(2J 2 



^) 4 



(65) 



<9r 



<9r 



is a constant of the motion, i.e., dE s /dr = 0. Notice that i? g is positive definite for any non 
trivial solution of the equations of motion. Therefore, it provides an absolute bound for each 
term in (165 p . Thus, if at any time r all the terms in ( 165]) are finite, then they will remain 
finite for any evolution that satisfies the condition of regularity for r = 0, demonstrating the 
absolute stability of the system under any finite perturbation for which all the terms in ( 165]) 
exist. Of course, a simple example is the case of an incoming wave of (essentially) compact 
support, for which we have already demonstrated the stability. 
We should further notice that the term, 



T,, 



(2J 2 + R 2 ) 2 d X i d X3 



R 4 



dr dr 



(66) 



r=R 



in the second line of (165]) allows for mild compensating divergences of dxi/dr, and dxs/dr, 
provided they are compatible with the existence of the integrals in ( I65p . We notice nev- 



ertheless that, taking into account (159]) . and 
forms, 



60]) . we may also write Td in the alternative 



AJ 4 (6R 2 J 2 + U 4 + R 4 ) 2 
Ro 10 {Ro 2 + J 2 ) 2 

4 J 4 (6R 2 J 2 + AJ 4 + R 4 y 
Ro W {Ro 2 + J 2 ) 2 



Xi (r, Rq) + 



J 2 (8R 2 J 2 



J 2 (4R 2 J 2 + 4 J 4 
Ro 3 (6R 2 J 2 + U 4 + R 4 ) 

+ 3R 4 + 4J 4 ^ 



£(t) (67) 



.Ro 3 {QR 2 J 2 + 4 J 4 + Ro 4 ) 



£(t)- X3 (r,Ro) 



and, since both Td and £(r) are bounded, we find that both Xi{ t iRq) an d Xz{ t iRq) must 
also remain bounded, and this implies that X2{t, r) is also bounded. On account of 052p . this 
result implies that j + (t,r) is also bounded to first order. We remark that 7_ is directly of 
second order and vanishes to first order. Therefore, the full metric remains bounded under 
arbitrary perturbations. 



XI. NUMERICAL SOLUTIONS OF THE INITIAL PLUS BOUNDARY 

CONDITIONS PROBLEM 

One can now check that the system (150]) . (155]) . (156]) . and (158]) . with appropriate replace- 
ments of xi by X3) can be used to set up and solve numerically the initial plus bound- 
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ary value problem where, for some arbitrarily chosen r = Tq we give arbitrary values to 
£0o), (dC(r)/dr)\ T=T0 , to Xifar) and (dxi(r,r)/dr)\ T =T0 , in the interval < r < i?o> 
and to X3( r 0; r ) an d (<9x 3 (t, r)/dr)\ T=T0 , in the interval _R < r, subject to the constraints 
dxi(T,0)/dr = 0, and those resulting from (I55|) . 

Although this procedure is well defined, and leads to a unique evolution for appropriate 
initial data, at least in its Cauchy domain, we must remark that the well-posedness of this 
type of problems, as regards sensitivity to small changes in the initial data, is still an open 
and interesting issue, outside the limits of the present research. The purpose of this exercise 
is mainly to display the evolution of some appropriately chosen initial data, in particular of 
£(t), in the region where we may expect the presence (or absence) of quasi normal ringing. 
As we shall see, this behaviour, at least in the region analyzed, is in complete agreement 
with our previous results. 

The integration procedure we have chosen uses a simple finite difference method for 
updating xi in < r < _R 0) an d for % 3 in R Q < r < r a , where r Q » R is some appropriately 
chosen outer boundary. The values of £ are updated using a simple leap frog scheme. Once 
these updates are carried out, we use the matching conditions (155]) to update Xii an d X3 
at r = Rq- We will be mainly interested in the behaviour of £ as a function of r. On this 
account we choose some appropriate value of r D and boundary condition for % 3 at r Q , but 
carry out the integration only from r = to r ~ ?"o/2, to ensure, by causality, that no signal 
coming in from r = r D has time to reach r = Rq and affect that behaviour. In particular, in 
all the examples below we set X3( r o,t) = 0, while we chose r Q = 50 for the first and second 
example, and r Q = 300 for the other two examples. 

The first explicit example of the results of this numerical integration for £(r) is given in 
Fig. 4 and Fig. 5. The initial data used was: 



and we set R = 0.25, J = 0.125, corresponding to x — 0.5, and r D = 50. In Fig. 4 
we display the evolution of £(r) as a function of r. The graph corresponds clearly to an 
evolution dominated by QNR, although a more detailed analysis indicates the presence of a 
small non oscillating background at large r. In Fig 5 we have a plot of ln(|£(r)|) (thick line) 
in the region dominated by the damped oscillation mode. The thin line curve corresponds 
to an approximate fit with the function ln(| cos(2.18r — 1.4) |) — 0.147r, although the actual 
plot in Fig. 5 has been displaced up by 0.5 for clarity. We can see again the very clear 
signal of a damped oscillation, with a frequency and fall off in very good agreement with 
those of the example of Fig. 2, although we have used a different procedure and form for 
the incoming pulse. The agreement might probably improve using a more elaborate rather 
than our simple minded implementation of the numerical procedure. In any case, it was 
not the purpose of the authors to optimize it but rather to show that even with a simple 
implementation we can solve the initial value problem related to the perturbation treatment 
of the dynamics of our system close to a static configuration, and obtain results in perfect 
agreement of those extracted from the mode expansion. 
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FIG. 4: £(r) as a function of r for initial data of the form (|68p . for i?o = 0-25 and J = 0.125, that 
is x = 0.5. 

As a second full numerical example we considered again the initial data given f!68p . but 
setting Rq = 0.25 and J = 0.25, that is x = 1.0. The results are shown in Fig. 6. In this 
case, in agreement with our analysis of the QNM, we find a very strongly damped oscillation, 
very much like that shown in Fig 3. for the same value of the parameters Rq and J, although 
the computational procedures applied in each case were completely different. We can also 
check that both the frequency and damping are in good agreement with the values given in 
Table I for x = 1.0 and positive real part for Q. 

As a final example, we present the numerically computed evolution of £(r) as a function 
of r for the cases Rq = 0.25 and J = 0.35, that is x = 1.4, given as the thin line curve in 
Fig. 7, and for R = 0.25 and J = 0.55, that is x = 2.2, also given in Fig. 7 as the thick line 
curve. In these examples the QNR is over damped and no trace of "ringing" is apparent in 
the time dependence of £(r). In particular, the last 2.2, corresponds to y = 4.84 in 

the notation of j6[ , where in what they call the first class solutions one has a pure imaginary 
QNM frequency, and, therefore, we would have only a purely exponentially decreasing QNR. 
These plots confirm also the stability of the corresponding static configurations. 



In this Section we consider the problem of the apparent discrepancy between the evo- 
lutions obtained using the formalism developed in |4j and those described here. We shall 
illustrate the reasons for this by considering the QNM toy model given by Price and Husain 
in [3] , since, as will be seen, it bear a certain formal resemblance to the shell model discussed 
here. The Price and Husain model consists of two torsion strings, the first attached at x = 
and extending to x = a, where it is attached to a second string that extends to x = oo. The 
strings have different characteristic constants, so that the respective equations of motion 
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FIG. 5: a = ln|f(r)| for initial data of the form ([BSD, for = 0.25 and J = 0.125 (thick line). 
The thin line is an approximate fit by ln(| cos(2.08r — 1.4) |) — 0.147r. Notice that the actual plot 
of this fit has been displaced up by 0.5 to avoid cluttering the figure. 



arc. 



and, 



d 2 ui d 2 u x 
+ 



dx" 1 



— k 



d 2 u 2 ^ d 2 u 2 







(69) 



(70) 



dt 2 dx 2 

where k > 1 is a constant. U\ = Ui(t,x), and u 2 = v,2(t,x), represent the torsion angle 
respectively in the intervals < x < a, and a < x < oo. The boundary and matching 
conditions are, 



Ml(t,0) 

ui(t,a) 
dui 



dx 





u 2 (t,a) 
du 2 



(71) 



dx 



The general solutions of (!6"9]) and (170]) . are, 

ui(t,x) = f(t + x)-f(t-x) 



(72) 



where /(?/) is, in principle, an arbitrary function of y, and we have included the boundary 
condition at x — 0, and, 

w 2 (t, s) = F(t + kx) + G(t - kx) (73) 
where and G(y) are also functions of y. 
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FIG. 6: £(t) as a function of r for initial data of the form ()68|) . for Rq = 0.25 and J = 0.25, that 
is x = 1.0. Notice the strong damping of the oscillations after the arrival of the incoming pulse. 




FIG. 7: £(t) as a function of r for initial data of the form (I68p . for i?o = 0-25 and J = 0.35, that 
is x = 1.4 (thin line), and -Ro — 0.25 and J = 0.55, that is x = 2.2 (thick line). We see that the 
quasi normal oscillations are absent or completely damped out, and we only have the response of 
the shell to the incoming pulse. 
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One can now check that, imposing the matching conditions, and assuming that all the 
functions vanish for large values of their arguments, the general solution functions f(y), 
F(y), and G(y) satisfy the relations, 

F{t + kx) = — [{k + I) fit + kx-ka + a)-(k- l)f(t + kx - ka - a)} 
G{t-kx) = — [(k-l)f(t-kx + ka + a)-(k + l)f(t-kx + ka-a)} (74) 

and, therefore, an arbitrary f{y) determines a complete solution of the problem. In partic- 
ular, if we choose f(y) of compact support, we will also get F and G of compact support. 
This type of solution describes then an incoming pulse of compact support that approaches 
the junction at x = a, excites a non vanishing u\ for a finite time, and gives rise to an out 
going pulse, also of compact support. There is, in general, no "ringing", as no quasi normal 
mode is excited. This is precisely the situation considered in jij, where we could fix the 
function £(r) arbitrarily, and, by choosing it of compact support we got also ip- and ip + of 
compact support, with no quasi normal ringing. 

But the Price and Husain model does display QNR. To see how this comes about we 
assume that the functions admit Fourier transforms of the form, 



f(y) = / e""f{u)du> 

J — oo 

/oo 
e iuiy F(uj)duj (75) 
-oo 

POO 

G(y) = / e iuiy G(cu)dcu 



(76) 

Replacing in the system (1741) we eventually find that we may write, 



"OO 



m = 2l e^ /T"yi k > 
Aoo e 2luja (l + k)-k + l 

G(y) = I" + ( 77) 

Kyj e 2 ^ a (1 + k) - k + 1 v ; 

These expressions provide the solution for the problem of finding the evolution assuming 
that we are given an incoming pulse F(t + kx) as the data for the problem, just as in the 
present shell problem. If F(u) has suitable analytic properties it may be possible to compute 
the integrals in (J77J) by extending the path of integration to the complex u plane. We notice 
here that the denominators in (1771) have simple poles at, 

7r i fk + 1\ /r7nX 

u = u n = n- + — In (78) 

a 2a \k — 1 J 

where n is an integer, including zero. Notice that all u n have the same positive imaginary 
part, and, therefore, when it is possible to close the path of integration by adding an infinite 
semicircle in the upper half plane of oo, each one of these poles will contribute a factor of the 
form exp(— ln((/c + l)/(k — l))y/ (2a)) times an oscillating factor depending on n and on the 
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explicit form of F(u). These are the characteristic features of the "ring down" associated 
with the quasinormal modes. In fact, if we go back to ( J74l) . and write it in the form, 

F(t) = — [(k + l)f(t — ka + a) — (k — l)/(t - ka - a)} 
G(t) = -r[(k- I) f{t + ka + a)-(k + I) fit + ka- a)} (79) 
and assume f(y) = exp(iuy) we get, 

F(t) = — [(k + 1) exp(iu;a) — (k — 1) exp(— iua)] exp(?u;(t — ka)) 

G(t) = — [(k — 1) exp(iua) — (k + 1) exp(— iua)} exp(?u;(t + ka)) (80) 

The condition Fit) = corresponds to a purely out going wave for x > a. This is achieved 
precisely for u of the form (178|) . Taken literally, these type of solution correspond to both 
f{t — x) and f{t + x) that decrease exponentially for t — > +oo, but increase exponentially as 
t — > — oo. Similarly, we have G{t — kx) that grows exponentially for large x. These solutions 
are therefore unphysical. What is then the relation between these solutions and the poles 
in f[77|) ? The crucial point is that the poles are effective only if the integration path can be 
closed on the upper half plane of u. In general, for k > 1, this requires y + Ika > 0, and 
therefore, the exponential terms are not present for y, and as a consequence t, less than a 
certain lower bound. 

But, why do we see a ring down in one case and not in the other? To understand what 
is happening here we go back to (!79|) . and assume that f(y) has a Fourier transform f(u). 
Then, replacing in (17^1) . we find, 

F(co) = f(u)e- iuka [(k + 1) e iua - e~ iua (k - 1)] (81) 

But this implies that for general f{oo), the transform F{ui) will contain a factor that precisely 
cancels the denominators in (177|) . and therefore there will be no ring down in the solution 
of the problem, unless, of course, f(u) itself contains the appropriate poles, as in (1771) . 
The general conclusion is then that both (17^1) and (JTTj) provide a complete solution of the 
problem and are, therefore, completely equivalent. Clearly, for the purpose of making the 
QNR apparent it is simpler and, in a sense, more "natural" to use (1771) and specify freely 
F(y), i.e., the incoming pulse shape, just as was done for the shell problem in the present 
work. 

XIII. FINAL COMMENTS 

In this paper we have considered again the perturbative analysis of the static configu- 
rations of the Apostolatos and Thorne shell model by modifying the formalism developed 
in jij. As a result we have been able to show the completeness of the mode expansion as 
regards the characteristic data problem for the perturbative dynamics of the shell. This, 
in turn, provides a simple and direct proof of the stability under bounded (symmetry pre- 
serving) perturbations of the general static configurations of the shell. We have also derived 
a set of coupled linear ordinary and partial differential equations that describe the general 
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perturbative evolution of the shell. This set of equation can be used to set up an initial 
value problem for the shell that can be solved numerically, but, more importantly, one can 
prove the existence of a positive definite constant of the motion, that implies the stability of 
the motion resulting for an arbitrary perturbation. In several examples considered we find 
perfect agreement between the full numerical evolution and that obtained through the nu- 
merical integration of the mode expansion, although they are completely different in detail. 
At first sight our results seem to be in contradiction with those of Kurita and Nakao M. We 
believe that there is no contradiction here, and that our computations and those in [6] are 
perfectly compatible. The problem appears because Kurita and Nakao conclude that the 
existence of certain complex zeros (or poles) in the complex uj plane automatically imply 
a direct effect on the evolution of the shell, but this is not necessarily the case. The usual 
arguments for relating these poles to the evolution relay rather heavily on the possibility 
of extending integrals on the real u axis to the complex u plane, in such a way that one 
effectively picks up the residues of those poles, but this may not always be possible, and, 
even, it might happen that one can disregard those poles, when they exist, by extending 
the integration path in the opposite direction. In fact, since the mode expansion involves 
Bessel functions in rather complex combinations, the matter of finding the extensions of 
appropriate function to the complex oj plane is a highly nontrivial undertaking. Going back 
to our own derivations, we notice that although we found QNM with a negative real part 
of u, we only find evidence of QNR for those with a positive real part for u. We remark 
once again that the mode expansion does not require explicit consideration of the complex u 
plane, as only the real u axis is involved. We believe that although the work j6[ is valuable 
as regards the finding of the QNM of the system, since these by themselves are unphysical 
(they diverge generally in some space or time direction), given the very complex nature of 
the functions and of their possible extensions to the complex plane, conclusions about the 
stability of the static configurations of the shell can only be reached by explicitly showing 
the relation between these modes and the general evolution of the shell. In fact, in view 
of our results concerning completeness of the mode expansion and its relation to stability, 
we can only conclude that all in principle unstable QNM must be suppressed, as well as 
possibly some of the stable ones. Of course, explicit confirmation of this conclusion would 
require performing and analyzing appropriate extensions into the complex plane, but that 
is completely outside the scope of the present research. 
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